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ABSTRACT 


This  paper  is  a  study  of  the  conditions  that  can  be  expected  for  the  in¬ 
tensity  of  light  reflected  and  scattered  from  the  earth.  The  albedo  of  the 
earth  is  defined  and  the  albedo  of  various  surfaces,  and  regions  of  the  earth 
are  discussed.  An  attempt  is  made  to  survey  typical  results  of  measurements 
of  Tiros  data.  The  albedo  of  the  earth  covered  by  a  Rayleigh  atmosphere  is 
presented  in  graphical  form  and  the  intensity  distribution  over  the  earth  is 
plotted  for  various  orientations  of  observer  viewed  spot  and  sun.  A  calcula¬ 
tion  of  expected  minimum  response  is  made  for  a  specific  photodiode  for  a 
certain  geometry  and  a  method  is  presented  to  enable  interested  persons  to 
easily  perform  similar  calculations  for  different  geometries. 
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SURVEY  ON  EARTH  ALBEDO 


Before  beginning  the  discussion,  the  oft-used  term  "albedo"  must  be 
clearly  understood.  The  albedo  of  an  object  is  simply  the  r<itio  of  the  emer¬ 
gent  flux  of  light  scattered  and  reflected  in  a  given  spectral  region  to  the 
incident  flux  of  light  in  the  same  spectral  region.  It  is  mathematically  com¬ 
putable  if  the  emergent  spectral  intensity  [in  units  of  (power)(unit  projected 
area)"^  (steradian)"^  (unit  wavelength  interval)  is  known.  Then 


ai =  ai  (xr  X2^ =  albedo  =  — 


»  X  o  r*  2  7T 

r7T  /2 

Ih  Jo 

I  (0 ,  cp)  cos  0  sin  0  d  0  d  cp 
J  o  X 

r  2 

Jxj  * 


F.  cos  Z  d  X 
A 


where 


I  is  the  emergent  spectral  intensity 

A. 

F  is  the  incident  spectral  intensity  =  incident 
K 

spectral  flux  [  (power ) (unit  area)  ^  (unit  wavelength  interval) 
divided  by  tt  for  a  source  like  the  sun. 

Z  is  the  solar  zenith  angle  shown  in  Fig.  L 
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The  albedo  just  defined  is,  strictly  speaking,  the  albedo  of  a  point.  The 
albedo  of  the  entire  object  would  have  to  be  obtained  by  averaging  MaM  over 
the  object. 

To  measure  the  albedo  of  a  point  on  the  earth's  surface  from  an  orbiting 
satellite  requires  quite  a  different  definition: 


a2  " 


cos  Z 


*  (X)d  X 


where 

F  is  the  flux  measured  by  the  sensor 

S  is  the  incident  spectral  flux 

$(\)  is  the  effective  spectral  response  of  the  sensor 

The  value  F  is  not  directly  measureable  because,  to  accurately  measure 
F,  sensors  would  have  to  be  placed  at  all  points  in  a  hemisphere  centered  on 
the  viewed  area  and  the  intensities  would  have  to  be  averaged  over  the  hemi¬ 
sphere.  Since  this  is  clearly  not  a  practical  idea, an  assumption  has  to  be 
made  as  to  the  scattering  function  of  the  area.  The  assumption  that  is  in¬ 
variably  made  is  that  the  surface  scatters  isotropically,  i.e.  ,  according  to 
Lambert’s  law.  It  is  then  a  simple  matter  to  produce  a  number  F  =  tt  I  since 
the  intensity  in  all  directions  is  assumed  constant. 

The  albedo  of  an  area  measured  in  this  way  is  dependent  (especially  for 
a  portion  of  the  earth  where  scattering  from  the  atmosphere  is  the  major 
contributor)  on  the  exact  shape  of  the  spectral  response  curve  of  the  sensor. 
To  see  this  assume  two  sensors  with  their  3  db  points  at  the  same  wavelength 
but  having  differently  shaped  spectral  response  curves  as  in  Fig.  2  below. 

If  the  relative  intensity  of  the  scattered  sunlight  as  function  of  wavelength  is 
as  shown  below  (Fig.  3),  sensor  (b)  will  record  a  much  higher  albedo  than 
sensor  (a). 
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X  (m>  X  t*j> 


Fig.  2 


Fig.  3 
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The  realization  that  the  measured  albedo  depends  on  the  exact  shape  of  the 
spectral  response  curve  as  well  as  the  bandwidth  has  never  been  encountered 
in  the  literature. 

Two  alternative  definitions  of  albedo  are  proposed  below,  although  they 
are  not  definitions  of  albedo  in  the  generally  accepted  sense  and  would  more 
properly  be  called  (for  want  of  a  better  name)  backscatter  coefficients. 

I,  (fi,cpf  Z) 

bn  (X)  - - - •=  spectral  backscatter  coefficient 

*  h> 

whe  re 

L  is  the  spectral  intensity  of  the  scattered  sunlight  at  a 

X 

given  wavelength 

H  is  the  spectral  intensity  of  sunlight  at  the  same  wavelength 
X 

Since  only  the  ratio  of  these  two  factors  is  important,  may  be  set  to 
1  for  theoretical  calculations.  The  spectral  flxix  seen  through  a  solid  "duo" 
from  a  point  above  the  surface  is  just  bjH^duo. 


b 


2 


,  cp,  Z)$  (X)dX 

— — F - =  sensor  backscatter  coefficient 

,  H  $  (X)dX 

«  A.  A. 


where  $  is  the  effective  spectral  response  of  the  sensor. 
The  effective  flux  flowing  into  "duo"  is  just 


b 


2 


Jx  Hx$  (X)duo  =  dS(0,  cp, 


Z) 


These  backscatter  coefficients  are  explicit  functions  of  3  variables  A,  cp, 
and  Z,  and  implicit  functions  of  the  type  of  area  viewed,  i.e.  ,  forests, 
deserts,  oceans,  etc.  ,  and  the  sensor  itself.  Hence,  they  would  be  very  difficult 
to  tabulate.  All  that  is  undertaken  in  the  latter  half  of  this  paper  is  to  obtain 
values  for  b^  and  b^  for  very  specific  minimum  intensity  cases  for  a  specific 
sensor . 

*  The  Fairchild  Planar  Photodiode  1N3734 
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The  following  is  a  survey  of  some  of  the  information  available  on  albedo 
of  the  earth  as  it  was  defined  in  this  paper. 

The  albedo  of  the  entire  earth  is  generally  computed  by  observing  earth- 
light  reflected  from  a  new  moon  or  by  satellite  measurements.  In  the  visible 
region  30#  -  40#  are  respectable  limits  for  the  albedo  of  the  earth  at  any  time 
with  a  time  average  value  of  36-39#  depending  on  the  investigator.  For  a 
cloudless  earth  the  albedo  is  16  ,  for  a  cloud -covered  earth  it  is  56#.  The 
albedo  is  about  50$  in  the  UV  region  and  28#  in  the  IR  region. 

The  albedo  of  various  cloud  types  for  the  entire  spectrum  is  shown  in 
Table  I  below.  ^ 


TABLE  I 

Albedo  of  Various  Cloud  Types. 

Cloud  Type  Albedo 

Very  dense  clouds  of  extensive  area  and  great  depth  78 

Dense  clouds,  quite  opaque  65-62 

Dense  clouds,  nearly  opaque  44 

Thin  clouds  36-40 

Str atocumulus ,  overcast  66-81 

Altostratus,  occasional  breaks  17-36 

Altostratus,  overcast  39-59 

Cirrostratus  and  altostratus  49-64 

Cir ros tratus ,  overcast  44-50 

Stratus,  680-1000  ft.  thick  78 


Data  for  several  types  of  terrain  are  given  below  in  Table  II. 
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TABLE  II 


Albedo  of  Various  Surfaces. 

Surface 
De  sert 

Fields,  various  types 

Forest,  green 

Grass,  various  conditions 

Ground,  bare 

Sand,  dry 

Sand,  wet 

Snow  or  ice 


Albedo  (per  cent) 


24-28 

3-25 

3-10 

14-37 

7-20 


18 

9 

46-86 


Water  (direct  sun  only)  Z  (degrees) 


0 

2.  0 

20 

2.  1 

40 

2.  5 

50 

3.  4 

60 

6 .  0 

70 

13.  4 

80 

34.  8 

85 

58.  4 

90 

100.  0 

where  Z  is  the  zenith  angle  of  the  sun. 

The  albedo  of  water  according  to  one  author  is  closely  approximated  by 

the  above  numbers,  calculated  from  Fresnel's  law  of  reflection,  even  when 

the  wind  is  10  mph.  Angstrom  states  that  in  ordinary  geophysical  problems 

the  data  in  Table  II  are  applicable.  ^ 

A  good  source  of  values  for  earth  albedo  is  the  results  of  radiation 

measurements  from  the  Tiros  meteorological  satellite  series.  The  results 

2 

gotten  by  Nordberg  and  Bandeen  are  given  below  as  an  example  of  what 
information  can  be  culled  from  Tiros  data. 
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The  Tiros  satellites  are  generally  equipped  with  two  types  of  radiometers. 
The  first  is  a  non- scanning  radiometer  with  a  broad  response  in  both  the 
visible  and  infrared  regions  and  a  rather  low  spatial  resolution.  Its  aperture 
is  approximately  55  degrees.  The  second  is  a  scanning  radiometer  which 
scans  as  the  satellite  spins.  It  is  of  medium  spatial  resolution,  has  a  5  de¬ 
gree  aperture,  and  responds  to  radiation  in  five  different  spectral  regions 
separated  by  filters.  Three  of  these  regions  are  in  the  infared  between  5.  0 
and  6.  7  microns,  between  7.  5  and  13.  5  microns,  and  between  7.  0  and  32.  0 
microns.  The  other  two  regions  lie  mainly  in  the  visible  portion  of  the  spec¬ 
trum,  that  is,  between  0.  2  and  7.  0  microns  and  between  0.  50  and  0.  75  mic¬ 
rons.  It  is  the  latter,  channel  5,  which  is  of  interest  here.  Note:  The 
satellite  is  flying  in  a  nearly  circular  orbit  with  an  approximate  altitude  of 
7  80  km. 

Shown  in  Fig.  4  through  Fig.  6  are  albedo  maps  made  with  the  channel  5 
sensor.  These  values  of  albedo,  as  all  other  albedo  values  derived  from  Tiros 
data,  assumes  an  isotropic  reflecting  surface. 

The  data  of  interest  for  minimum  intensity  considerations  is  to  the  left 
of  the  sub- satellite  path  where  specular  reflection  from  the  water  does  not 
occur.  Minimum  albedo  for  this  area,  with  the  sun  near  zenith,  is  4"\ 
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Figure  4.  Tiros  Albedo  Percentages 

Albedo  percentages  measured  by  channel  5  (approximately 
0.  5  -  0.75 fj.  )  of  the  medium  resolution  radiometer  between 
1439  and  1445  GMT,  Orbit  117,  Z0  July  1961.  The  skies 
are  mostly  clear,  except  for  some  minor  scattered  clouds. 
Averages  of  these  data  are  given  for  the  areas  enclosed  by 
dashed  lines.  The  heavy  dashed  line  is  the  sub-satellite 
path  with  points  each  minute. 
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The  next  albedo  map  is  for  the  eastern  United  States  under  mostly  cloudy 
skie  s. 
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Figure  5.  Albedo  Percentages  Measured  by  Channel  5. 

Notice  the  low  value  of  8^  in  the  absence  of  clouds. 
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The  last  albedo  map  shows  the  North  African  Desert  under  clear  skies. 


|  3-63-~ 39?  I 


Figure  6.  Albedo  Percentages  Measured  by  Channel  5. 

The  albedo  values  presented  in  the  above  maps  do  not  take  into  account  the 

fact  that  the  earth  does  not  scatter  isotropically.  To  get  a  quantitative  idea 

as  to  the  extent  to  which  anisotropy  occurs,  the  following  computation  with 

3 

Tiros  III  data  was  performed. 
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A  selected  area  of  the  earth.  Fig.  7,  is  viewed  by  the  radiometer  twice 
during  the  orbit  with  the  scans  3  to  4  minutes  apart.  The  seitellites  have 
moved  far  enough  in  the  time  for  the  scattering  angle  of  the  scattered  and 
reflected  radiation  to  have  changed  quite  significantly.  The  two  maps  above 
show  the  same  area.  Notice  the  difference  in  the  appearance  of  the  two  maps, 
a  difference  which  is  not  acknowledged  in  the  previous  albedo  maps.  It  should 
be  noted  that  the  smaller  scattering  angles  in  this  study  are  associated  with 
the  larger  nadir  angles.  There  is  a  large  difference  in  nadir  angle  between 
the  two  sets  of  data  along  X.  However,  when  this  difference  is  reduced,  as 
along  Y,  the  fluxes  for  the  smaller  scattering  angles  remain  considerably 
higher.  Since  the  solar  zenith  angle  varies  in  an  identical  manner  along  both 
lines,  the  difference  in  flux  can  be  least  qualitatively  be  attributed  to  the  dif¬ 
ference  in  scattering  angle.  The  variation  of  channel  5  flux  along  X  and  Y  for 
the  two  different  viewing  geometry  is  shown  in  Fig.  8. 

The  data  of  both  solar  channels  suggest  that  the  scattering  angle  effect  is 
largest  in  areas  of  least  cloudiness.  This  suggests  the  possibility  that  spec¬ 
ular  reflection  from  the  ocean  surface  could  have  been  responsible  for  the 
larger  flux  differences  in  the  areas  of  relatively  low  flux.  Along  Y,  this  effect 
could  have  entered  since  this  line  lies  nearly  in  the  plane  of  the  sun.  Since 
both  X  and  Y  show  the  same  relationship  between  flux  difference  maximas 
and  regions  of  relatively  low  flux,  specular  reflection  alone  could  not  have 
been  the  cause.  Increased  atmospheric  scattering  for  the  smaller  scattering 
angles,  regardless  of  the  viewed  surface,  is  probably  the  important  factor. 

In  the  report  just  considered,  the  data  sample  included  significant  differ¬ 
ences  in  the  radiometric  nadir  angle  associated  with  the  separate  scans  of  the 
overlapping  data.  Therefore,  the  increase  in  forward  scatter  could  have  been 
effected  by  the  variation  in  the  viewed  volume.  However,  even  as  in  a  sub¬ 
sequent  publication,  when  the  nadir  angle  difference  of  the  two  scans  is  re¬ 
duced  to  zero,  the  measured  albedo  remain  higher  for  smaller  scattering 
4 

angles  (Fig.  9). 

Notice  in  the  figure  that  in  general,  effective  blackbody  temperature 

measured  with  the  channel  2  sensor  (approximately  8.  0-12.  0  microns)  decreases 

^scattering  angle  180°  -(angle  between  sun,  viewed  spot  and  satellite  sensor) 
**nadir  angle  ~  angle  between  local  vertical  at  the  sensor  and  the  direction  of  the 
viewed  spot. 
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Figure  7.  Relative  Patterns  of  Channel  5  Flux  for  Area  with  Overlapping 
Scanlines  from  Orbit  5,  Tiros  III 

(a)  Initial  Alternating  Open  Mode  Data 

(b)  Complemental  Single  Open  Mode  Data 
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Figure  8. 


Variation  of  Channel  5  Flux  for  Radiometer  Scans  Towards  and 


Away  from  the  Sun  and  the  Variation  of  Their  Difference  Along 
X  and  Y. 
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Figure  9.  Comparison  of  Data  from  the  Wall  and  Floor  Sensors  Along  the 
Isopleth  of  Zero  Nadir  Angle  Difference 

(a)  Channel  2  Effective  Temperatures 

(b)  Channel  5  Albedos  and  Albedo  Differences 
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as  the  albedo  increases.  This  is  due  to  the  fact  that  the  tops  of  clouds  are 
much  colder  than  the  less  reflective  underlying  surface. 

It  would  seem  that  a  careful  study  of  raw  Tiros  data  could  yield  information 
on  the  minimum  light  intensities  and  reflected  and  scattered  from  the  earth. 
This  path  of  investigation  was  not  chosen,  however,  because  of  two  consider¬ 
ations.  The  spectral  regions  of  the  Tiros  sensors  and  the  sensor  under  con¬ 
sideration  are  different  enough  to  make  the  Tiros  data  almost  inapplicable  for 
this  study;  also  it  has  been  noticed  that  many  authors  tend  to  disregard  Tiros- 
derived  albedos  whenever  they  fall  below  the  albedo  that  can  be  expected  from 
scattering  by  a  Rayleigh  atmosphere.  As  an  example  one  author  found  the 
albedo  of  a  point  over  ocean  to  be  6#.  For  the  same  geometry  and  for  approxi¬ 
mately  the  same  spectral  regions  the  albedo  of  a  standard  atmosphere  with 
Rayleigh  scattering  and  zero  surface  reflectivity  deduced  by  computations  by 
K.  Li.  Coulson  is  8^.  It  is  stated  that  the  difference  is  dependent  on  the  de¬ 
gradation  of  the  instrumental  response,  although  a  proper  evaluation  should 
include  the  effects  on  the  Rayleigh  albedo  of  molecular  absorption  and  the 
scattering  and  absorption  by  aerosols. 

RAYLEIGH  SCATTERING 

The  problem  of  atmospheric  scattering  is  not  an  easy  one  and  has  been 

5 

tackled  only  recently  by  Chandrasekhar  in  his  book  on  Radiative  Transfer. 

The  equations  which  he  introduces  have  been  tabulated  by  K.  L.  Coulson  in 
his  various  publications  on  scattering  by  a  Rayleigh  atmosphere.  Some  of 
his  results  and  the  results  of  others  working  along  similar  lines  are  discussed 
in  the  succeeding  pages. 
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The  values  that  follow  have  been  computed  for  parallel  unpolarized 
radiation  of  net  flux  tt  per  unit  area  incident  on  the  top  of  an  idealized  plane 
parallel  atmosphere.  The  atmosphere  is  assumed  to  be  uniform  and  infinite 
in  the  horizontal,  non-absorbing,  composed  of  scattering  centers  which  are 
small  compared  with  the  wavelength  of  the  light,  and  of  finite  optical  thichness. 

For  the  sun  at  a  given  zenith  angle  0^  the  total  outward  flux  from  the  top 
of  the  atmosphere  in  a  unit  frequency  interval  and  unit  time  is  given  by 


D 

v 


2tt  7t/2 

J  «. 

o  o 


I  (P,  cp)  sin  0  cos  0  d  0  d  (p 


0,  ,  and  are  defined  by  Fig.  1  with  Z  replaced  by  6^.  I  (0,  cp)  is  a 

tabulated  function. 

The  spectral  albedo  or  the  monochromatic  reflectance  is 
D 

R  = - - - 

v  7 rp 

o 

where 

u  =  cos  0 
o  o 


Below,  R  ^  is  tabulated  for  various  wavelengths  and  zenith  angles. 

Table  III  is  represented  graphically  in  Fig.  10.  The  values  indicate  the 

reflectance  increases  for  decreasing  wavelength  at  all  sun  elevations  due  to 

5 

the  inverse  fourth  power  dependence  of  Rayleigh  scattering. 
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TABLE  III 


Reflectance  of  the  Atmospheric  Model  at  Various  Wavelengths  for  Representative 
Values  of  Sun  Elevation  and  Surface  Reflectivity.  (Lambertian  reflector)',:t> 


A 

- o — 

X3150A 

4000 

5000 

6000 

7000 

8000 

0 

0.  02 

0. 

754 

0. 

638 

0. 

56  5 

0. 

503 

0. 

430 

0. 

325 

0 

0.  10 

0. 

705 

0. 

552 

0. 

396 

0. 

247 

0. 

154 

0. 

094 

0 

0.  40 

0. 

550 

0. 

304 

0. 

152 

0. 

079 

0. 

044 

0. 

024 

0 

0.  80 

0. 

400 

0. 

181 

0. 

085 

0. 

039 

0. 

022 

0. 

013 

0 

1.  00 

0. 

345 

0. 

148 

0. 

068 

0. 

033 

0. 

017 

0. 

009 

0. 

25 

0.  02 

0. 

791 

0. 

712 

0. 

66  5 

0. 

625 

0. 

574 

0. 

494 

0. 

25 

0.  10 

0. 

752 

0. 

643 

0. 

536 

0. 

427 

0. 

362 

0. 

319 

0. 

25 

0.  40 

0. 

623 

0. 

445 

0. 

345 

0. 

298 

0. 

276 

0. 

285 

0. 

25 

0.  80 

0. 

492 

0. 

347 

0. 

293 

0. 

271 

0. 

260 

0. 

257 

0. 

25 

1.  00 

0. 

450 

0. 

322 

0. 

280 

0. 

264 

0. 

256 

0. 

255 

0. 

80 

0.  02 

0. 

922 

0. 

910 

0. 

904 

0. 

896 

0. 

885 

0. 

86  5 

0. 

80 

0.  10 

0. 

918 

0. 

892 

0. 

866 

0. 

842 

0. 

826 

0. 

815 

0. 

80 

0.  40 

0. 

860 

0. 

830 

0. 

812 

0. 

805 

0. 

803 

0. 

800 

0. 

80 

0.  80 

0. 

811 

0. 

798 

0. 

798 

0. 

798 

0. 

798 

0. 

799 

0. 

80 

1.  00 

0. 

794 

0. 

792 

0. 

794 

0. 

796 

0. 

798 

0. 

799 

Note  a  more  complete  table  by  Coulson  is  given  in  Appendix  I- 
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[3-63-3925 


WAVELENGTH  (A) 


Figure  10.  Total  relative  monochromatic  flux  emerging  from  the  top  of  the 

atmosphere  as  a  function  of  wavelength,  for  various  sun  elevations 
and  three  values  of  surface  reflectivity:  A  =  0  (solid  curves), 

A  =  0.  2  5  (dashed  curves)  and  A  =  0.  80  (dotted  curve  s).  Curve  s  are 
labeled  with  values  of  ju  .  A  curve  (X)  for  primary  Rayleigh, 
scattering  and  no  attenuation  is  included  for  comparison. 
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This  graph,  however,  is  only  relative  and  assumes  all  the  incident  flux 
to  be  the  same  at  all  wavelengths.  If  we  assume  a  source  with  the  same 
variation  with  wavelength  as  the  sun  the  graph  of  Fig.  11  is  obtained. 


3-63-3926 


Figure  11.  Absolute  monochromatic  flux  from  the  top  of  the  model  atmosphere 
as  a  function  of  wavelength,  for  zero  surface  reflectivity  and  five 
different  sun  elevations  (ju  =  1.  0,  0.  40,  0.  20,  0.  10,  0.  02;  A  =  0) 


Using  the  graphs  of  Fig.  11  the  albedo  for  a  plane  parallel  earth  as  a  function 
of  sun  elevation  for  the  spectral  region  3200  A  ^  X  ^  8000  A  can  be  obtained. 

The  results  obtained  by  Coulson^  are  plotted  in  Fig.  12  and  tabulated  in  Table  IV. 
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TABLE  IV 


o  o 

Values  and  Comparisons  of  the  Total  Absolute  Flux  fqr  3200  A  \  <  8000  A 
at  Six  Different  Sun  Elevations.  (Units  of  flux:  W/m  ) 


Sun 

elevation 

1.1° 

5.  7° 

11.5° 

23.  6° 

53.  1° 

90.  0° 

F 

o 

15.  4 

77.  0 

154.  0 

307.  9 

615.  8 

769.  9 

(J)o 

8.  0 

25.  0 

34.  3 

44.  0 

48.  9 

51.  8 

(J>25 

9.  7 

37.  1 

62.  4 

103.  6 

179.  9 

214.  8 

(J>80 

13.  7 

65.  9 

128.  3 

249.  6 

488.  6 

609.  6 

<">o 

0.  519 

0.  325 

0.  223 

0.  143 

0.  079 

0.  067 

(«')25 

0.  630 

0.  482 

0.  405 

0.  336 

0.  292 

0.  279 

^80 

0.  890 

0.  856 

0.  833 

0.  811 

0.  793 

0. 792 

<F0g>0 

(F0g>25 

(F0k)80 

7.  4 

52.  0 

119.  7 

263.  9 

566 .  9 

718.  1 

7.  6 

53.  2 

122.  1 

272.  4 

581.  2 

740.  1 

8.  5 

55.  5 

128.  5 

291.  5 

636  0 

801.  5 

(*) 


A 


Incident  flux  on  horizontal  surface  at  top  of  atmosphere. 

Emergent  flux  through  horizontal  surface  at  top  of  atmosphere. 
(Subscript  A  represents  value  of  surface  reflectivity:  0.025  or 
0.  80). 

(J)A 

— -  -  Albedo:  ratio  of  emergent  to  incident  flux. 

0 


(F 


0g)A 


Flux  incident  on  horizontal  surface  at  ground. 
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Figure  12.  Approximated  albedo  for  the  entire  solar  spectrum  (solid  curve) 
and  computed  albedo  for  the  spectral  range  1200  A  >  X  ^  8000  A 
(dashed  curve)  as  a  function  of  sun  elevation  (A  =  0). 

Coulson  also  did  a  calculation  of  the  planetary  albedo  of  the  spherical 
earth  due  to  atmospheric  scattering.  The  disk  of  the  earth  was  divided  into 
concentric  plane  parallel  rings  and  the  sun  elevation  for  each  ring  was  con¬ 
sidered  constant.  A  value  of  7.6f  was  obtained  for  the  entire  spectrum.  This 

o 

value  is  reduced  by  ozone  absorption  in  the  region  near  and  below  3000  A. 

The  corrected  albedo  for  this  more  realistic  case  is  given  as  6.  9$. 

7 

Some  additional  graphs  computed  by  Coulson  in  a  later  publication  for  a 
plane  parallel  atmosphere  with  a  Lambertian  surface  at  the  bottom  are  pre¬ 
sented  in  Figs.  13  -17.  Figures  13-15  are  a  graphical  representation  of  the 
values  tabulated  in  Appendix  I.  It  can  be  easily  seen  that  the  atmosphere  it¬ 
self  scatters  back  a  large  portion  of  the  short-wavelength  radiation  especially 
for  the  larger  solar  zenith  angles,  and  that  at  longer  wavelengths  and  smaller 
solar  zenith  angles  the  effect  of  the  reflectance  of  the  underlying  surface 
becomes  predominant.  Figures  16  -17  show  the  strong  dependence  of  system 
reflectance  for  large  solar  zenith  angles  particularly  pronounced  for  small 
values  of  A  (for  the  entire  solar  spectrum  (Fig.  16)  and  for  3  spectral  regions 
(Fig.  17)). 
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MONOCHROMATIC  ALBEDO  OF  PLANE  PARALLEL  ATMOSPHERE 


Fig.  13  Monochromatic  reflectance  of  model  atmosphere -plus  - 
surface  system  as  a  function  of  surface  reflectivity  for  optical 
thickness  t  =  1.  00  for  five  different  solar  zenith  angles  j  =  1.  00, 
corresponding  to  wavelength  X  =  3120  A  at  sea  level. 
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MONOCHROMATIC  ALBEDO  OF  PLANE  PARALLEL  ATMOSPHERE 


Fig.  14  Monochromatic  reflectance  of  model  atmosphere-plus- 
surface  system  as  a  function  of  surface  reflectivity  for  optical 
thickness  t  =  0.  10  for  five  different  solar  zenith  angles  t  =  0.  10, 
corresponding  to  wavelength  \  =  5460  A. 
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Fig.  15  Monochromatic  reflectance  of  model  atmosphere -plus  - 
surface  system  as  a  function  of  surface  reflectivity  for  optical 
thickness  t  =  0.  001  for  five  different  solar  zenith  angles  t  =  0.  001, 
corresponding  to  wavelength  t  =  17,  000X. 
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SYSTEM  REFLECTANCE 


3-63-393 1~1 


Fig.  16  Reflectance  of  model  atmosphe  re -plus -  surface  system 
as  a  function  of  cosine  of  solar  zenith  angle  (j,  for  spectral  interval 
3100  A  \  25,  000  K  for  diffe  rent  values  of  surface  reflectivity. 
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Fig.  17  Reflectance  of  model  atmosphe re -plus- surface  system  as 
a  function  of  cosing  of  solar  zenith  angle  ju  ,  with  spectrum  divided  into 
ultraviolet  (3100  X  <  X  <  4000  X),  visible  (X000  X),  and  infrared 
(7000  A  <  X  <  25,  000  A) 
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The  graphs  and  tables  that  have  been  presented  so  far,  for  the  most  part, 
have  indicated  the  total  amount  of  light  reflected  from  the  plane -parallel 
atmosphere.  Relatively  little  information  has  been  presented  about  the  ex¬ 
pected  minimum  light  conditions  at  various  configurations  of  observer,  viewed 
spot  and  source.  The  mater>.l  on  the  succeeding  pages  will  provide  information 
on  the  expected  intensities  and  nbackscatter  coefficients11  of  the  reflected  and 
scattered  light. 

Some  of  the  results  are  presented  in  a  form  which  makes  it  necessary  to 
discuss  briefly  the  concept  of  optical  thickness,  t.  t  is  defined  as 

z 


T 


'  O 

n(z) 

o 


ct.a(z)  dz 


where 

n  =  number  of  scattering  centers/unit  volume 
=  scattering  cross-section 

is  usually  taken  as  the  height  of  the  atmosphere 

Then  the  intensity  of  light  on  the  ground  is 


I =  I  e“T 
o 

t  is  related  to  \  by  the  following  Table  (Table  V)  and  Fig.  18. 


TABLE  V 


1.  00 

.  50 

.  25 

.  15 

.  10 

.  05 

.  02 

3120 

2715 

4365 

4950 

5460 

6440 

8090 

All  graphs  in  Figs.  19  and  20  take  into  account  the  effect  of  primary  and 
multiple  scattering  (unless  otherwise  indicated)  and  increased  path  lengths 
for  oblique  rays.  The  relative  intensities  indicated  on  the  graphs  are 
relative  to  an  incident  intensity  of  1.  0. 
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OPTICAL  THICKNESS  (t) 


Fig.  18  Optical  thickness  of  the  earth's  atmosphere  as  a  function 
of  wavelength. 
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90  70  50  30  10  10  30  50  70  90 


NADIR  ANGLE  (deg) 

<p  ~  -  0  (p  ~  p  -  7T 

0  0 


Figure  19.  Relative  intensity  for  three  different  values  of  optical  thickness 
as  a  function  of  direction  in  the  sun's  vertical.  Dashed  curves 
are  primary  scattering  only;  sc^lid  curves  include  higher  order 
scattering  (0  =  66.  3°,  A  =  0). 
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3-63-3935 


Figure  20.  Relative  intensity  as  a  function  of  direction  in  the  sun's  vertical, 

for  moderate  optical  thickness,  sun  elevation  53°  and  three  values 
of  surface  reflectivity  (t  =  0.15,  0^  =  36.  9°). 

The  plot  in  Fig.  20  shows  that  for  moderate  optical  thickness  the  surface 
reflected  light  may  constitute  the  major  part  of  the  emergent  light.  For 
values  of  A  approaching  1  the  surface  begins  to  look  like  a  Lambertian  surface. 
For  higher  optical  thickness  the  reflectivity  is  relatively  less  dominant  but 
very  important,  even  for  optical  thickness  near  1.  0. 

g 

In  a  recent  paper  by  William  Snoddy  the  case  of  the  plane-parallel  atmos¬ 
phere  was  extended  for  a  spherical  earth.  The  portion  of  the  earth  visible 
at  a  particular  altitude  was  divided  into  400  plane  parallel  sections  and  the 
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intensity  computed  for  each  section.  The  irradiation  incident  on  a  differential 
element  at  a  given  altitude  was  computed  and  compared  to  the  irradiation 
expected  if  the  earth  were  a  diffuse  surface  to  obtain  an  ‘’effective  albedo". 

Note  that  this  definition  is  different  from  the  other  definitions  of  albedo  pre¬ 
sented  earlier  and  could  have  extensive  use  with  an  earth  sensor  of  very  wide 
aperture.  The  geometry  of  the  situation  is  presented  in  Figs.  21-22. 

Towards  the  limb  of  the  earth,  incident  light  travels  through  a  smaller 
air  mass  than  if  the  atmosphere  were  plane-parallel.  The  air  mass  through 
which  an  incident  beam  will  travel  has  been  calculated  by  Bemporad  and  the 
air  mass  is  given  in  Fig.  23  together  with  the  air  mass  for  a  plane  parallel 
atmosphere.  This  graph  is  used  to  correct  the  data  obtained  for  a  plane 
parallel  atmosphere  for  large  angles  of  incidence. 

Figures  24  -  56  are  plots  of  intensity  over  the  planetary  sphere  with 
isophotes  drawn  through  the  points  of  equal  intensity.  The  sun  moves  to  the 
right  along  the  horizontal  diameter  of  the  circles.  The  differential  element 
is  above  the  center  of  the  circles. 

Figures  24  -  30  illustrate  the  variation  as  a  function  of  0  with  H  =  10^  km, 
o  o 

T  =  1.  0  (X  =  3120  A)  and  A  =  0.  Figures  31  -  37  show  the  same  variation  with 

3 

0  but  have  H  =  10  km.  Variation  of  intensity  with  altitude  is  illustrated  in 
o 

Figs.  38  -  42.  Maps  showing  the  effect  of  the  reflectance  of  the  underlying 
surface  are  given  for  an  altitude  of  a  million  kilometers  (Figs.  43  -  45)  and 
a  thousand  kilometers  (Figs.  46  -  48). 

Perhaps  the  most  interesting  maps  from  the  point  of  view  of  using  a  light 
sensor  with  a  broad  rather  than  a  line  bandwidth  are  the  maps  showing  vari¬ 
ation  in  intensity  with  optical  thickness  and  hence  wavelength.  Examples  are 
shown  for  T  =  1.  0,  0.  5,  0.  25,  0.  1,  0q  =  30°;  A  =  0  and  H  =  10*'  km  (Figs.  49  - 
52)  and  10*  km  (Figs.  53  -  56). 

These  relative  values  of  intensity  are  what  has  been  defined  previously  as 
the  spectral  backscatter  coefficient,  "b^".  The  quantity  bjH^  da)  was  then 
integrated  over  the  sphere  and  the  relative  flux  intensities  plotted  in  Figs. 

57  -  58.  As  the  reflectance  of  the  underlying  surface  is  increased,  the  maximum 
flux  intensity  shifts  away  from  the  ultraviolet  towards  the  maximum  of  the  sun. 
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The  spectral  effective  albedo  of  the  earth  is  shown  in  Fig.  59  for  H  =  10  *  km 

and  three  values  of  A.  It  should  be  realized  that  this  graph  is  valid  only  for 

the  given  altitude.  In  Figs.  60  -  6l  the  effective  albedo  for  the  spectral  region 

o 

between  3050  and  8050  A  (graphs  1)  and  for  the  entire  spectrum  (graph  2) 

were  calculated  as  a  function  of  A  and  H  for  two  values  of  0  . 

o 
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EARTH 

Fig.  21  Section  on  planet  surface 
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Fig.  22  Element  above  planet 
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AIR  MASS  (0J/AIR  MASS  (0  =  0) 
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Fig.  24  Map  of  isophotes  for  H  =  10  km,  0  =  0 
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Fig.  28  Map  of  isophotes  for  H 


10  km, 
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Fig.  29  Map  c£  isophotes  for  H 


10  km,  0  =13  5° 
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Fig.  30  Map  of  isophotes  for  H 


10^  km, 
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Fig.  36  Map  of  isophotes  for  H  =  10  km,  0q  =  90° 
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Fig.  37  Map  of  Isophotes  for  H  =  10  km,  0q  =  95° 
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Fig.  39  Map  of  isophotes  for  H  =  10  km 
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Fig.  40  Map  of  isophotes  for  H  =  10  km 
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Fig.  41  Map  of  isophotes  for  H  =  10  km 
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Fig.  42 


Map  of  isophotes  for  H  =  10  km 
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Fig.  43  Map  of  isophotes  for  H  =  10  km,  A  =  0 
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Fig.  44  Map  of  isophotes  for  H  =  10  km,  A  -  0.  25 
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Fig.  45  Map  of  isophotes  for  H  =  10  km,  A  =  0.  80 
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Fig.  51  Map  of  isophotes  for  H  =  10  km,  j  =  0.  25 
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Fig.  52  Map  of  isophotes  for  H  =  10  km,  t  =  0.  10 
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SPECTRAL  EFFECTIVE  ALBEDO 
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Fig.  59  Spectral  effective  albedo  for  A  =  0,  0.  2  5  and  0.  80; 

0=0°  and  85°. 
o 


71 


EFFECTIVE  ALBEDO 


Fig.  60  Effective  albedo  as  a  function  of  altitude  for  Pq  =  0° 
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EFFECTIVE  ALBEDO 
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Continuing  the  type  of  calculations,  the  results  of  which  were  presented 
above,  the  intensity  distribution  seen  by  a  sensor  at  synchronous  altitude  is 
presented  below  for  specific  cases.  The  numbers  have  been  arrived  at  by 
using  an  extremely  useful  reference.  "Tables  Related  to  Radiation  Emerging 
from  a  Planetary  Atmosphere  with  Rayleigh  Scattering"  by  Coulson,  Dave 
and  Sekera. 

In  what  follows  it  has  been  assumed  that  the  sun,  the  satellite  earth  sensor, 
and  the  center  of  the  earth  all  lie  on  a  straight  line.  The  symbols  used  are 
defined  by  Figs.  21-  22.  The  sensor  sweeps  from  a  spot  directly  below  to 
the  edge  of  the  earth.  Corrections  have  not  been  made  for  the  curvature  of 
the  earth's  atmosphere  because  corrections  are  only  necessary  for  very 
oblique  angles  of  incidence  (sun  zenith  angles  over  82°)  and  in  such  cases 
the  scattered  intensity  is  large  (for  this  geometry). 

Figure  62  illustrates  the  intensity  (relative  to  an  incident  intensity  of 
1.  0)  of  scattered  light  as  a  function  of  angle  6  with  wavelength  as  a  parameter. 
The  same  information  is  presented  in  Fig.  63,  but  here  the  intensity  is  a 
function  of  the  sensor  nadir  angle.  This  makes  the  sharp  increase  in  intensity 
near  the  limb  of  the  earth  much  more  apparent.  The  increased  backscatter 
near  the  limb  is  the  reason  why  a  sharp  halo  of  violet  rings  the  earth  in  the  color 
Gemini  IV  photographs. 

Although  the  relative  intensity  of  scattered  light  increases  as  wavelength 
decreases,  the  intensity  of  sunlight  decrease  s  as  X  approaches  UV  wavelengths. 
The  graphs  of  Fig.  64  have  been  multiplied  by  the  spectral  flux  of  the  sun  to 
obtain  the  relative  intensity  of  scattered  sunlight  at  various  wavelengths. 

The  data  contained  in  these  graphs  has  been  multiplied  by  the  sensor 
response  at  two  values  of  6  and  displayed  as  a  function  of  \  in  Fig.  65.  Note 
that  the  graph  on  the  right  for  6  =  78.  45°  is  two  times  larger  than  the  graph 
for  6-0.  The  graphs  have  been  put  in  proper  perspective  in  Fig.  66  where 
graph  (2)  is  for  6  =  78.45°  graph  (3)  is  for  6=0°  and  graph  (1)  is  for  the 
intensity  distribution  over  wavelength  for  a  diffusely  reflecting  atmosphere¬ 
less  earth  for  6  =  0. 
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INTENSITY  (intensity  of  sun  radiation  is  1,000  at  oil  X's) 


Fig.  62  Relative  intensity  of  scattered  light  at  synchronous 
altitudes,  as  a  function  of  the  central  angle  of  a  great  circle 
($  =  180°,  ot  =  180°,  0o  =  0°). 
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INTENSITY  (intensity  of  sun  radiation  is  1,000  at  all  X's) 
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X  INTENSITY,  tw/cm /u/*teradian) 


Fig.  64  Relative  intensity  of  scattered  sunlight  at  various  wave¬ 
lengths  as  a  function  of  the  central  angle  of  a  great  circle 
($  =  180°,  O'  =  180°,  eo  =  0°). 
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x  INTENSITY,  x  SENSOR  RESPONSE  (8*0)  (x  10 


WAVELENGTH  (p) 


Fig.  65  Light  intensity  seen  by  sensor  after  adjustment  for 
sensor  response. 
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Fig.  66 
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RELATIVE  ENERGY 
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The  area  under  the  curves  is  proportional  to  the  flux  seen  by  the  sensor 
and  the  current  produced.  The  earth  has  been  assumed  flat  and  the  intensity 
constant  in  the  region  viewed.  The  ratio  of  the  area  of  graph  (2)  or  (3)  to  the 
area  of  graph  (1)  is  the  sensor  backscatter  coefficient  T,t>^!t. 

To  predict  the  current  produced  by  the  sensor  a  very  simple  experiment 
can  be  performed.  Measure  the  current  produced  by  the  sensor  looking 
straight  down  at  the  light  reflected  from  a  diffusely  reflecting  flat  source, 
illuminated  by  one  solar  constant  with  the  same  wavelength  distribution  of 
intensity  as  the  sun.  If  this  current  is  I  and  if  the  response  of  the  sensor 
is  linear,  the  current  output,  I,  in  outer  space  for  a  particular  configuration 
is : 


I  =  b9I 
2  o 

b^  for  6  =  0.  0  is  .  0201 
b^  for  6  =  78.  45  is  .  204 

The  calculation  was  done  as  follows:  The  geometry  of  the  sun,  the  earth,  and 
the  satellites,  is  drawn.  A  value  of  y  is  chosen  and  the  corresponding  value 
of  6  is  gotten  from  Appendix  II.  Then  y  +  6  =  cos  V*  .  The  angles  cos 
and  cp  are  gotten  from  the  geometry.  The  relative  intensity  can  then  be  found 
in  Table  VI  of  Coulson’s  tables  for  all  r  and  hence  X«  The  value  of  intensity 
at  a  particular  X  is  multiplied  by  the  flux  from  the  sun  at  the  same  wavelength 
in  watts/cm  /micron.  (The  relative  spectral  distribution  of  the  sun  is  given 
in  Fig.  6  7.  The  sun's  output  has  been  smoothed  out  so  that  undue  weight  is  not 
placed  on  dips  in  the  curve.  This  curve  must  be  multiplied  by  .  220  to  make 
the  curve  absolute.  )  This  value  is  then  multiplied  by  the  relative  sensor  re¬ 
sponse.  The  data  points  are  then  plotted  versus  wavelength  in  microns  and 
the  area  under  the  curve  is  computed.  The  ratio  of  this  area  to  the  area 
under  curve  (1)  in  Fig.  66  is  the  sensor  backscatter  coefficient  ,,b>n.  This 
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_4 

area  is  189.  70  boxes  or  474.  25  x  10  sq.  units  (i.  e.  each  box  is 

_4 

.  05  x  .  005  =  25  x  10  sq.  units).  Alternately,  if  the  curve  (1)  is  multiplied 
by  the  relative  sensor  response  and  the  area  is  found  the  answer  is  the  same 
but  the  work  is  reduced. 

Time  did  not  permit  additional  calculations  of  b  7  for  different  satellite 

orientations.  It  is  felt  that  the  minimum  possible  value  of  b  >  for  any  geometry 

is  not  much  lower  than  .  02.  The  calculation,  however,  does  not  take  into  account 

the  effects  of  water  vapor  and  molecular  absorption  of  radiation  at  large  values 

of  wavelength.  Ozone  absorption  plays  a  minor  role  since  it  becomes  important 
o 

only  below  3000  A.  Scattering  by  dust,  haze  and  aerosols  would  serve  to 
increase  the  intensity  of  scattered  light  beyond  that  of  a  Rayleigh  atmosphere. 
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APPENDIX  I 

The  Monochromatic  Reflectance  Rx  of  the  Atmosphere-Surface  System  for 
Various  Values  of  the  Optical  Thickness  r  and  Surface  Reflectivity,  A. 
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APPENDIX  I 


The  Monochromatic  Reflectance  Rv  of  the  Atmosphere -Surface  System  for 
Various  Values  of  the  Optical  Thickness  t  and  Surface  Reflectivity,  A. 
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APPENDIX  II 

Angles  y  and  6  Defined  by  Fig.  21  for  the  Radius  of  the  Earth  =  6371  km 
and  the  Radius  of  the  Satellite  Orbit  =  42,  16 0  km.  Angles  are  Given  in 
Degrees. 
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